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ABSTRACT 



We are considered the propagation of the fundamental and second har- 
monic solitary waves under the slowly varying envelope pulses approximation 
in a quadratic nonlinear medium that characterized by negative refraction 
index at the frequency of fundamental wave and by positive refractive index 
at the second harmonic frequency. We find the some solutions of the evo- 
lution equations which described the steady state simulton, coupled second 



O ■ harmonic and fundamental frequency waves propagation in this materials. 



PACS: 42.65. Tg 
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^ : 1 Introduction 



In the past few years, new developments in structured electromagnetic ma- 
terials have given rise to negative refractive index materials which have both 
negative dielectric permittivity and negative magnetic permeability in some 
frequency ranges. These materials are often referred to as left-handed mate- 
rials (LHM) or materials with negative refraction (NRIM). The properties 
of such materials were analyzed theoretically by Mandelschtam and Veselago 
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many years ago d I2]- It is very important that the wave vector is anti- 
parallel with the Poynting vector inside LHM PP-[7]. The occurrences of the 
LHM in microwave region have been demonstrated experimentally in the [S]- 
jllj . The experimental evidence for the existence of NRIM in optical region 
was found [12]- [1^] also. The results of investigation and discussion of the 
applications of negative refractive medium represent in reviews |17| ll8j. 

Unusual features of the NRIM manifest itself in the passing the interface 
between this material and the positive refractive index one. Also, the re- 
fractive index of the same medium can be positive in one frequency region 
and to be negative in another frequency region. Hence the features of NRIM 
would be made itself evident in interaction of wave packets having carrier 
frequencies from different spectral regions, where refractive indexes have op- 
posite signs pHl 1^ . Second harmonic generation (SHG) is first example 
where negative refraction leads to distinctly different pictures of spatial dis- 
tributions of the interacting waves intensity [23 I2n]- Here the quadratic 
nonlinear NRIM acts similar to distributed Bragg reflector, i.e., Bragg grat- 
ings. Namely, the harmonic wave propagates towards fundamental (pump) 
wave. In nonlinear Bragg grating under some conditions the connected pair 
of the incident and reflected waves develops. It is gap soliton propagating 
in nonlinear grating [2Zj. We can expect that in quadratic nonlinear NRIM 
the second harmonic and fundamental wave pulses are able to produce the 
complex steady state solitary wave that is similar to quadratic soliton [^ 1^ 
or simulton [2H]- 

In this paper the propagation of the complex steady state of the funda- 
mental wave and second harmonic wave considered. It is assumed that refrac- 
tive index for the frequency of fundamental wave is negative, and refractive 
index at the harmonic wave frequency is positive. The particular solutions 
of the equations describing the evolution of these waves are founded. They 
are analogues of the solitons [211 122] and cnoidal waves [23 [^ in quadratic 
nonlinear medium. Under condition that value of the phase mismatch is very 
large the propagation interacting waves approximately governs by nonlinear 
Schrodinger equation [211 1321 • 



2 Constituent equations 

Consider the propagation of the pulses of parametrically coupled waves, i.e., 
fundamental wave (or pump wave) and second harmonic one. The type 
I phase-matching condition is taking place PHI 1201 • In this case the slowly 
varying envelopes of pump and harmonic waves evolve according to following 
system of equations can be written in following form 



^[ki^ + -Tu]E,-^—E^ = l^lL^p^^(cui)exp(-zA;i4l) 



dz vidtj 2 dt"^ c^ki 

where ui = ujp,U2 = oug = 2ujp, k^ = {Ljj/cYe{ujj)fi{ujj), kj is sign of square 
root of k^, and 



Pnl{uji) = Xal^i; ^2, -uJi)E2El eyiY>[iz{k2 - fci)], 
Pnl{'^2) = X2i^2;^i,^i)E'^exp[iz{ki + ki)]. 

where D12 are dispersion coefficients at pump and signal wave frequencies |2H 

122] • Here Ei is slowly varying envelope of the pump (fundamental frequency) 
wave, and E2 is the slowly varying envelope of the second harmonic wave. 

As we have in equations ()ll2p ratio {uj j / c) ^J ^{uj j) / e{(jj j) , the coefficients 
in right part of equations ()1|2|) are not depending on signs of the dielectric 
permittivity and magnetic permeability. Nevertheless direction of the wave 
propagation is defined by signs kj. 

It should be pointed that we consider the collinear propagation and so 
the collinear phase mismatch. The condition Ak = means that vector 
k2 of pump wave is collinear and equal to vector 2ki. The configuration of 
these vectors defines the energy fiow directions. If the nonlinear medium 
is characterized by positive refraction index, then vectors ki^2 and Pointing 
vectors have same directions. Opposite figure take place if the nonlinear 
medium is characterized by negative refraction index. In this case waves with 
frequencies from frequency region of negative refraction index are directed 
according to zero phase mismatch condition, but the Pointing vectors are 
contra-directed [231221 (see Fig.l). 



We will consider the SHG under assumption that frequency of pump wave 
belong the frequency region of negative refraction index, while the frequency 
of the signal wave lies in frequency region positive refraction index. 

Let us introduce the normalized variables 

C = z/L, T = {t + z/vi)/tp, El = Aiod, E2 = ^2062 exp{izAk), (3) 



^10 = c{2niUiLx2ii^i)y^Veiuj^yj4uJi), (4) 

A20 = c{27ruj2Lx2iuJ2)y^\/^i^2)/K^2), (5) 

Xi = x^^n^i;2t^i7-^i)> X2 = x^^\'^^i;^i,^i), (6) 

/^l,3 = /^(^l,3), £^1,3 = £^(^1,3), (7) 

where A A; = k2 — 2ki. Then SHG is described by the following normalized 
system of the equations 

■d a d^ * n fQ\ 

'dc'' + 2 9^"^ - '''' = ° (^^ 

where 



-1/..-1 , ..-i\ ^ n o ( -Ai \ h{^i)e{2uji 




L<,->r' + ^f ') > 0. »=|tVI r^ii^^ij^i^r' 



a^agnDi, ,8 = ZJj/IAI, L = ti/\Di\. 



pi 



Parameter 5 takes account of the walk-off effect for pump and harmonic 
pulses that is due to the difference of the group velocities directions for 
interaction waves. This parameter can not be zero, as in the case of positive 
refractive index medium. 

3 Integral of motion - the Manley-Rowe re- 
lation 

By using the equations (J8I9|) the following expression can be obtained 
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From these equations it follows that 

(3 d I ^de2 de*2\ a d f ^dei de\ 

If we take into account the boundary condition 

W.2? -^ |eio,2or = const, 



or 



at \t\ -^ ±oo, then 



lei.r 



e2p — T^kiP 1 '^T = const. 



In the case of the SHG in continuum wave that results in the Manley- 
Rowe relation jSEj 

|e2p kiP = const. 

When SGH takes place in positive refractive index medium the minus in 
this expression is replaced to sign plus. It corresponds to the propagation of 
the both interacting waves in one direction. 



4 Real form of the equations 

It is suitable to introduce the real variables for the interacting waves 

ei = aexp(i(y9]^), 62 = 6exp(iy92)- (10) 

Substitution of these expressions into ()8|9|) leads to 



da a f da du), d'^Lo. 



dC 2\ dr dr dr"^ 

TrT + 5- ^2- — ^ + h ^ 



dC, dr 2 \ dr dr dr"^ 

(Jj I ~T ~ (Jj 



dC 2 \dT^ dr dr 



a6sin$, 


(11) 


1 2 . ^ 
-a sm$, 


(12) 


— a6cos$, 


(13) 


1 2 ^ 
-a cos$, 

2 


(14) 



d(^ dr J 2 ydr"^ dr dr 

where ^ = (f2 ~ 2(^9^. The case of continuum wave has been considered in 
fIF\, EE] 5 where principal difference between spatial distribution of the pump 
and harmonic fields in the NRI medium with respect to positive refraction 
index medium was pointed. 

5 Steady state solutions 

Let us consider solutions without chirp, i.e., ones with d'^ip^ 2/^'^'^ = 0- 
Hence, we have dipi 2/^'^ — ^1,2- The mismatch condition $ = 0, which 
is assumed additionally, leads to ^^2 = 2r2i = 2Q. Furthermore, we assume 
that dipi .2/^C = '*i,2-, thus we suppose the phase functions evolve as 

ip^ = k( + Qt, <P2 = 2kC + 2Qt. (15) 

The amplitude equations (fTT|) and (fT^ take the following form 

ac+""97 = 0, (16) 

|+(^-2,«n)f ^ 0, (17) 

The steady state waves for both frequencies must propagate as single one. 
Thus, from ^T^ fTTjl we get the amplitudes as function of single variable 

y = r-C/V 

a = a{T-(/V), b = b{T-(/V) 

with 
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This expression defines the parameter "instant frequency" Q = S/{a + 2(3). 
The phase equations (jl3|) and (J14j) take the following form 

— - - {2aK + n^)a = 2aab, (18) 

^ + 1(A + 2k + 2SQ- 2f]Q^) b = \a\ (19) 

OT^ (j ^ P 

Let us suppose that b = fa. From (fTS|) and (fT^ we get two equations for 
one function a = a{y) of one variable y = t — (^/V : 

-^~(2aK + n^)a = 2afa\ (20) 

^ + l{A + 2K + 26n-2pQ')a = i-a^ (21) 

dy^ (3 ^ (3f 

These equations would define single function a, if the coefficients before equal 
orders of a are same, i.e., p = a/ 2(3 > and 

- (2//5) {A + 2K + 26n - 2l3n^) = {2aK + fi^). (22) 

It should be remarked that condition a^ = 1 was used. This expression 
defines the parameter (effective wave number) k 

?,(3Vl? - A5VL - 2A 
2a{2a + (3) 

From equation for function a{y) (J20|) we can obtain 

{da/dyf = {2aK + n^)a^ + (4(T//3)a^ + const. (23) 

If 0" = +1 , then dispersion parameter of second harmonic (3 must be 
positive too. Else, if cr = —1, then parameter (3 must be negative. 

Let be (T = —1. In this case we have / = (1/2|/5|)^/^, and the equation 
can be rewritten as 

{da/dyf = {n^ - 2K)a'^ - {4f/3)a^ + const. (24) 

Let us introduce p = (f2^ — 2k) and w{y) = {4:f/3p)a{y). Consider the 
boundary condition as following: a — > ao, da/dy — i> at r ^ ±oo. When 



Oq = 0, we will get the solitary wave on zero background. The constant 
of integrating is defined now by following expression const = —p{wq — Wq), 
where wq = (4//3p)ao . From equation (j^ the equation for w{y) follows 

{dw/dyf = p{w - Wo) [A^ ~{w- 7o)'] , (25) 

where 

A^ = (1 - wo){l + 3wo)/A, 7o = (1 - ^o)/2. 

If we introduce new variable u through the formula w = Wq + u'"^ , then 
the equation (|25j) transforms in to 

(du/dy)' = (p/4) [A V - (1 - j,uY] , (26) 

where 7^^ = (1 — 3wo)/2 . Parameter A^ is positive at < Wq < 1. 

5.1 Bright soliton solution 

The solution of this kind corresponds with choosing wq = 0, da/dy — >■ at 
r — > ±cxo. In this case the parameters in (j26p are A^ = 1/4 and 7;^ = 1/2. 
The equation (j26|) transforms into following one 

{du/dy)' = {p/4){u'-l), (27) 

Put u = cosh -(9, then from this equation we get {d'd/dy) = {p/4:) > 0. 
Hence, solution of (|7fj) is 



u 



{y)=cosh[p'/\y-yo)/2], (28) 



where yo is constant of integrating. Steady state solution of equations (fTTl 
inj) in the case of a = — 1 corresponding with bright soliton is 



(3p/4)V^ 
cosh'^[^{y-yo)/2] 

(3p/4) 
cosh^[^{y-yo)/2] 

If p = —\p\ < 0, equation (jTTj) has periodical solution 



<y) = „„';[':;:'^ '"',,„. . (29) 

Kv) = ...... 'r^V w„, - (30) 



u{y) = cos vb|(y-Z/o)/2 



which leads to singular solution of the system (fTTl -IT ^ . 

If a = +1 , then / = (1/2/3)^/2 ^ ^j^e equation ^^ takes the form 

{da/dyf = (fi^ + 2K)a^ + (4//3)a^ + const. (31) 

By introducing new function w{y) = — (4//3p)a(y) from (pTT|) the equation 
for w{y) follows 



(dw/dy) = p{w — Wo) A — {w — 7q) 
where p = (fi^ + 2k) , wq = — (4//3p)ao , and 

A' = {l-wo){l + 3wo)/^, 7o = (1 - ii}o)/2. 
Hence, this case reduced to the one considered above. 

5.2 Cnoidal wave solution 

The solution of this kind corresponds with choosing < wq < 1. If < 
Wq < 1/3 then we have A > 7^^ > 0. Equation (j26|) can be written as 

(du/dy)' = ip/A) [{A + 7,)«2 _ 1] [(^ _ ^^y + ^ (33) 



Introduce the new variable (p according to expression cosh(/) = Ma/A + 7^. 
It leads to 

{d<i)/dyf = (p/4) ( A + 7i) [1 + m^ cosh^ 0] (33) 

where m^ = (A — 7i)(A + 71)^"'^. Substitution cj) = iip into this equation 
leads to 

{di)/dyf = -(p/4)(A + 7i)(l+m^) [l - m^ sin^ 7/-] . 

It is equation for the Jacobi elliptic functions. Hence, we can write the 
solution by following formula 

sin?/' = sn(z2;,m), (34) 

where z = a/p(A + 7i)(l + w?){y — yo)/2, m^ = m^(l + rri^)^^. By using 
the properties of the trigonometric functions and Jacobi elliptic functions, in 



particular, sn{iz, k) = isc{z, k'),k''^ + k"^ = 1 , and sn^(2;, k) + cn^(2;, A;) = 1, 

we can write _ 

sn( z,m') 



sinhi 



cn[ z,m 



Now, we have 



u\y) 



(A + 7i)cn2(2;, m') 
Thus, the solution of the equation (^^j) represents by formula 

w{y) = Wo + {A + 7i)cn^(z, m'). 

By using definitions of the parameter fh and m^ one can find that 



(35) 



(36) 



m 



2 V A 



1/2 



2pA 

A + 7i 



1/2 



(y - yo) 



At interval < Wq < 1/3 module of Jacobi elliptic functions is not equal to 
unit, hence only cnoidal wave exists at these values of the pedestal amplitude 

Wq. 

The solution of the equation (j^Uj) under condition that 1/3 < Wq < 1 can 
be found in a similar way. If 1/3 < i^o < 2/3 this solution represents by 
formula 

w{y) = wo + {A~ |7j)cn2(zi,mi), 

where 

V_A^ '/' 

2 



(37) 



mi 



2 V A 



1/2 



Zl 



(y-vo)- 



If 2/3 < Wo < 1 we can find 



f \ , /A , I IN cn^ (2:2, "^2) 
w{y) = wo + {A+ |7il)r^7 



sn''[Z2,m2J 



(38) 



where 



m2 



2A 



nV2 



LA+I71I 



Z2 = l[\p\{A+\j,\)f'{y-yo). 



This solution describes the cnoidal wave, however amplitude of this wave 
periodically attains the infinite value. Hence, this solution is not acceptable 
for considered there physical situation. If wq > 1 equation (f^ has only 
singular (unlimited) solutions. 
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5.3 Dark soliton solution 

The quadratic dark soliton has been found in |331 ^^ . The similar solution 
can be found in the NRI medium. Let consider the equation ()25|). Substitu- 
tion in w = Wo + q results in the following equation 

{dq/dyf=pq[A'-iq-^,y]. (39) 

By following [33] let us find the solution in the form q = A/ cosh ay. If 
substitute that in (jHUj) . and equating the coefficients of various powers of 
cosh ay in this equation to zero, we obtain the following system of relations 

p(A2 - 7^) = 0, 2p7iv4 = Aa'^A, A^ = AaA. (40) 

From these relations it follows that a^ = p'j-^/2, A = 4a, and A^ — 7^ = 0. 
By using the definition of A^ and 7^ we can find that equation A^ = 7^ will 
be held at wq = and wo = 2/3. In the past case 7^^ = —1/2. From relation 
a^ = P7i/2 it follows that p must be negative, and a = ±(|p|/4)-^/^. Hence, 
we have A = ±2{\p\y/^. Thus, 



±Vbl 



cosh'[^\y/2]' 
The solution of equations (jllI - IT^ in this case is 

a{y) = (3|p|/4)v^ (2/3 ± 2v^sech^[/H?//2]) (41) 

b{y) = (3b|/4) (2/3 ± 2v^sech2[v^y/2]) (42) 

These expressions describe the quadratic dark soliton propagating in nonlin- 
ear medium with negative refractive index. 

5.4 Two-hump soliton 

Up to this point the solutions of equations (jTHj) and (fT^ have been found 
under constraint b = fa. It is known [35J that for positive refractive index 
medium there are solutions where this constraint is absent. These solutions 
were named by simultons in ^H]. We can show that the similar solutions 
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exist in NRIM too. Follow "SF we assume that solution of the (|18|) and (|19|) 
has the form 

a{y) = A tanh ay sechay, (43) 

b{y) = Bsecli^ay, (44) 

where A, B and a are unknown parameters. Substitution of this anzats in 
flTHj) and (Uni), and grouping the terms with the equal powers of cosh ay and 
tanhaysechay in these equations to zero, provides the following system of 
algebraic relations 





a^A-piA = 0, 


(45) 




2aAB + 6a^A = 0, 


(46) 




A^lS-Qa^A = 0, 


(47) 




-A^p + Aa^B = p2, 


(48) 


where {2aK + fi^) = 


= pi and (2//?) (A + 2fi; + 26n - 2/3^2) , 

a^ = p^> 0, A^ = Gpa^B, 
B = p2/2a^ B = -3aV. 


= P2- It leads to 



From these expression we obtain A"^ = — 18/5a^cr . It means that if cr = — 1 
then P > . Else, we can chouse cr = 1, then /? < and p2 = —\p2\ < . 
Let be 0" = — 1. In this case we can found 

a^=pi, B = 3pi, A = 3^/2(3pi, (49) 

and relation 

P2 = Qpl- (50) 

From (|MH) it follows 

( A + 2k + 25n - 2/3^2) = 3/3 {2k - J7^) ^ (51) 



The parameter f2 has been defined as fi = 5/ {213— 1), hence relation ()5Tj) 
defines parameter through phase mismatch A. From ()5H) it follows 

(±) _ 1 + 6/30^ ± y/l + 3/3 (A + 2(3/3 - 1)(]2 + 9/3(]4) 
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Thus, the two-hump sohton has the following form 



a{y) = 3pi V 2/5 tanh {^/jhy) sec h (y/ply) , (53) 

b{y) = Spisech.'^i^y). (54) 

For example, let us assume that /3 = 1, then A = —5^, 2k^'^^ = 5^ ± 1/3 and 
p[ = ±1/3. Only p[ = 1/3 is acceptable, hence, we get the two-hump 
solution of the system (fT8|l and (fT9|l 

a{y) = V2tanh(y/V3)sech(y/V3) , (55) 

b{y) = sech^fy/Vs). (56) 

For the case of the positive refractive index medium there is the solu- 
tion describing the perimetrically connected waves of the pump and second 
harmonic. Analogue of that exists in the NRI medium too. 

6 Conclusion 

We considered the steady state propagation of the coupled pair of waves, 
i.e., the pump wave at the frequency Uq and wave at the second harmonic 
frequency, in the quadratic nonlinear NRIM. In this case the group velocities 
of pump and harmonic pulses are directed in opposite directions j2Sl 1211 
123 I2n]- We are foud that under some conditions the self trapping of the 
interacting wave packets takes place. Two frequency wave packet (simulton) 
is generated from initial pulses alike to quadratic soliton of positive refractive 
medium (PRIM). 

To do comparison of the solitons in quadratic NRIM with solitons in 
PRIM it is suitable to represent the appropriat expressions for second case. 
The system of equations describing the quadratic solitons propagation in PRI 
medium can be written as 

• ^ ^ ^^ , * n ^t;7^ 
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where S = Ltp^{v^ ^ ~ "^2 ^) = tp\Di\ "^(f i ^ — ^2 ^)- In NRIM parameter S is 
always positive. Hence, there is no solutions corresponding with solitons of 
PRIM which have 5 = 

Velocities of the bright solitons in NRIM and PRIM are differed. For 
NRIM all steady state pulses move with velocity 



1 1 /^^ 2/5\ 1 fDi 2D 



Vs a + 2p\v2 vij Di + 2D2\V2 vi 
Whereas in PRI medium bright soliton has the velocity 



(59) 



1 1 ni-^)= ^' ^ fggl-£l|. (60) 



Vs 2(3 - a \vi f 2 y 2D3 - Di \ vi f 2 

In the plane {6, A) the existence regions of the steady state wave are 
differed also. For example, if one put dispersion parameters as follows a = 
/? = — 1 in the case of NRI medium from ()29IIH()|1 the expressions for envelopes 
of bright soliton result in 

a{y) = , , (61) 



6 cosh' 



l^{3A-S')/2{y-m)/3] 



,„) . JBZ£1 . (62) 



6cosh2[y/(3A-52)/2(y - yo)/3] 



This solution exists under condition A > 5 /3 . The solution of the equation 
(|Frj) and (J3H)) yields expression for quadratic soliton in PRIM 

a{y) = , , (63) 



2 cosh^ 



V(A + 5^)/2(y - 2/0)] 



KV) m^ • (64) 



2 cosh^ 



[^iA + 6')/2iy - yo)] ' 



In this case the solution exists under condition A + 5^ > . Thus the NRIM 
solitons and PRIM solitons exist in different regions of parameters plane 
(5, A). Stability conditions we plane to consider in future investigation. It is 
possible that the stability regions will be different too. 
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